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Abstract 

We investigate the interplay between topological charge and the 
spectrum of the fermion matrix in lattice-QED2 using analytic meth- 
ods and Monte Carlo simulations with dynamical fermions. A new 
theorem on the spectral decomposition of the fermion matrix estab- 
lishes that its real eigenvalues (and corresponding eigenvectors) play 
a role similar to the zero eigenvalues (zero modes) of the Dirac oper- 
ator in continuous background fields. Using numerical techniques we 
concentrate on studying the real part of the spectrum. These results 
provide new insights into the behaviour of physical quantities as a 
function of the topological charge. In particular we discuss fermion 
determinant, effective action and pseudoscalar densities. 
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1 Introduction 



Topological concepts play a prominent role in our understanding of contin- 
uum gauge theories. An example of particular importance is the idea that for 
the true vacuum functional of QCD (6'-vacuum) gauge field configurations of 
all topological charges have to be taken into account in the path integral [||]. 
Such statements however are difficult to formulate in a mathematically pre- 
cise way. The topological charge is a concept for classical, i.e. differentiable 
fields, since only those fields can be uniquely classified with respect to an 
integer topological charge. On the other hand differentiable configurations 
are of measure zero in the path integral for the gauge fields (see e.g. [0). 

When formulating gauge theories on the lattice the notion of a classical 
gauge field configuration gets lost. However, various definitions of a topo- 
logical charge for lattice gauge fields have been proposed (see e.g. the review 
[§] and also |0] for an overview on more recent developments). Using such a 
prescription one can in principle assign a topological charge to every lattice 
configuration (up to so-called exceptional configurations which in the con- 
tinuum limit do not contribute to the path integral). In a certain sense this 
makes the lattice approach very powerful for pursuing topological concepts 
since the path integral can indeed be decomposed into topological sectors 
whereas in the continuum formulation topological arguments are restricted 
to (non-contributing) classical configurations. 

The working hypothesis, that every gauge field configuration contribut- 
ing to the lattice path integral in the continuum limit can be assigned a 
topological charge, allows to decompose the lattice path integral into topo- 
logical sectors. As discussed above, in the continuum this is only possible 
in a formal way. The topological decomposition of the lattice path integral 
is a conceptionally sound procedure and allows to address many interesting 
physical questions in a mathematically meaningful setting. 

For classical gauge field configurations the Atiyah Singer Index Theorem 
provides a powerful tool connecting the topological charge with the zero 

modes of the Dirac operator. However in the continuum path integral its 

value is reduced to formal or semi-classical arguments, since the classical 

configurations are of measure zero as remarked above. 

On the lattice the situation is different. On one hand when approaching 

the continuum limit, every gauge field configuration contributing in the con- 
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tinuum limit can be assigned a topological charge. On the other hand we 
have no such thing as the Atiyah Singer Index Theorem on the lattice. How- 
ever, in the literature several investigations - |TT| attempting to establish 
the Index Theorem on the lattice in a probabilistic sense can be found. If 
such a goal is achieved, together with the fact that the lattice path integral 
can be decomposed into topological sectors, this would allow to gain better 
insight into the behaviour of fermionic observables using the connection be- 
tween topology and the spectrum of the fermion matrix. 



In this article we study QED2 on the lattice to demonstrate that a tight 
connection between topology and spectrum of the fermion matrix can be 
found and provides a powerful mathematical tool. Establishing this con- 
nection is a two step procedure. In a first step we use analytic arguments 
to understand general features of the spectrum of the fermion matrix. In 
particular we prove a new theorem on the chiral properties of the eigenvec- 
tors of the fermion matrix, which makes clear that only the real eigenvalues 
play a special role similar to zero eigenvalues (and zero modes) of the con- 
tinuum Dirac operator. These analytic results can be generalized to the 
4-dimensional Wilson-Dirac operator |jl2[. 

Once the real eigenvalues are identified as being the trace of the topolog- 
ical charge in the spectrum, in a second step we concentrate on analyzing 
their properties using numerical methods. Tracing only the real eigenvalues 
is a much simpler task than investigating all eigenvalues of the fermion ma- 
trix. For sufficiently large (3 we establish results on size and distribution of 
the real eigenvalues and show that their number is related to the topolog- 



ical charge computed with the geometric definition |T^. In a final section 
we make use of these results and use them to discuss the behaviour of the 
fermion determinant, the effective action and the pseudoscalar densities. 



The model under consideration (QED2) resembles many features of 4- 
dimensional gauge theories such as QCD. In particular U(l) gauge theory in 
2 dimensions allows for classical topologically nontrivial configurations ('vor- 
tices') which play the role of the instantons of 4D Yang-Mills theory. QED2 
has an anomaly which is related to the mass generation of the pseudoscalar 
singlet as in QCD. Thus the lattice version of this low dimensional model 
is an interesting candidate when exploring the role of topology in lattice 
models. Moreover for the continuum model there exists a limit which can 
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be solved analytically. When setting the fermion masses to zero, the model 
reduces to the Schwinger model |14| for two species of fermions. This allows 
to check the lattice simulations and to gain intuition from the analytic result. 
For the case of massive fermions semiclassical results and expansions in 
the bare mass are available |T6|. 



The article is organized as follows: Sec. 2 contains a short discussion of 
the model, the details of the simulation as well as definitions and theorems 
for the topological charge in the continuum and on the lattice. In Sec. 3 
we develop the analytic results for the spectral decomposition of the fermion 
matrix and in Sec. 4 we analyze the spectrum using numerical techniques. 
Finally in Sec. 5 we discuss the applications outlined above. The article ends 
with a discussion (Sec. 6). 



2 Setting 

This section prepares the ground for the program outlined in the introduc- 
tion. We formulate the lattice model, present technical details of the simula- 
tion and discuss the definitions and results for the topological charge in the 
continuum as well as on the lattice. 



2.1 QED2 on the Lattice 

We work on a two dimensional lattice A with volume L?. Lattice sites are 
denoted as x = (xi, X2) with Xi = 1, 2, L. The lattice spacing is set equal 
to 1. The gauge fields are group elements Uu{x) G U(l) assigned to the links 
between nearest neighbours x,x + and their action is given by 

^5 = /^E [1 - Re[/p(x)] , (1) 
where the plaquette element is defined as 



Up{x) = Uiix) U2ix + 1) Ui{x + 2) U2{x) . (2) 

The gauge fields obey periodic boundary conditions f/y(L + l, 0:2) = Uy{l,X2)i 
f/j/(xi, L + 1) = Up{xi, 1). The fermion action is the bilinear form (in matrix 
notation) 
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where the (imphcit) summation is over lattice points and spinor indices - 
the above mentioned two flavors come about through squaring the fermion 
matrix (see below). We write the fermion matrix as 

M = 1-kQ , (3) 

where the hopping matrix Q is defined as (spinor indices still suppressed) 

Q{x,y) 

Q^{x, y) = (1 + 7^)C/i.(x - i')A^_o,y + (1 - 'yu)U,y{x)A^+0^y . (4) 



i/=l,2 



j3 is related to the bare coupling constant e via (3 = 1/e^ and the hopping 
parameter k, to the bare mass m through k = (2m + 4)~^ The matrices 
are chosen as the Pauh matrices = a^-, v = 1,2,3. The fermions obey 
mixed boundary conditions, i.e. periodic in 1-direction and antiperiodic in 
2-dircction. This is taken into account by the mixed periodic Kronecker delta 
defined as 



x-l,y 



x+l,y 



x-2,y 



''x+2,y 



X2,y2 



■'X2,y2 



■'XI, yi 



^xi,l^yi,2 + •••• 

— ^X2,l^i/2,L + ^X2,2^y2,l + •••• 



XI, yi 



'X2 



lSyi,L + Sxi,2Syi,l + .... Sxi,LSyi,L-l 
Sxi,L-lSyi,L + Sxi,L^yi,l 
^y2,l + •••• ^X2,LSy2,L-l 
l^y2,2 + ■■■■ 5x2,L-l5y2,L - Sx2,LSy2,l 



(5) 



Obviously A 



y—i/,x 



^x+i',y 



We are rather explicit here, since the ana- 
lytic results for the spectrum of the fermion matrix in a fixed background 
configuration depend on the boundary conditions and also change whether 
L is even or odd. The latter dependence can be discussed conveniently us- 
ing Ax±j}^y. It has to be remarked, that when actually integrating over all 
gauge configurations the boundary conditions for the fermions are irrelevant 
for this model, since gauge group U(l) contains also —1. For a given gauge 
field configuration there is a partner differing just in the boundary links by a 
factor —1. The operators resulting from the Grassmann integral over fermion 
thus are averaged over these configurations. 
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2.2 Technical details of the simulation 



We use lattices of size with L between 4 and 16. In the Monte Carlo 
simulation the fermions were considered with the hybrid Monte Carlo method 
r^. In order to ensure positivity of the fermion determinant measure we 



had to introduce a second species of fermions. We stress that we squared the 
fermion determinant only for generating Monte Carlo configurations, whereas 
below we analyze the spectrum and other properties of M as defined in (|^) 
and (^) (not M^). In the HMC method we used trajectories with 10 steps 
and a step size adjusted such that the acceptance rate in the Monte Carlo 
step was 0.8 in the average. 

For the inversion of the fermion determinant even-odd preconditioning 
was used, together with the BiCG75 (actually BiCGds) algorithm [Q. This 
proved to be efficient for most configurations. Whenever we found slow con- 
vergence (so-called "unstable configurations" , also often called exceptional, in 
contradistinction to the notation for configurations without definite topologi- 
cal charge to be discussed later) we switched the inverter and used BiCGStab 
(biconjugate gradient stabilized, cf. e.g. without preconditioning, which 
seemed to be the most reliable (although not most efficient) method. 

In general we do not observe any abundance of unstable configurations 
close to Kcrit but actually only very few (0(0.1%)), which then were success- 
fully inverted by BiCGStab. In fact, we even can work at values above the 
presumed Kcrit without noticeable problems. 

In the computations below, we simulate the model at the values of k 
given in Table 1. These numbers have been chosen from preliminary results 
on KcritiP), that have been determined with help of restoration of PCAC 
along the lines of EO]. (The final results and their size dependence are dis- 



cussed in . It turns out that the k- values of Table 1 are slightly (typically 
0.001) below the critical values determined in |2T|.) These values are consis- 
tent with the behaviour of other observables like the chiral susceptibility. We 
used these k- values for all lattice sizes. As will be discussed below, our results 
are only weakly dependent on the precise value of the hopping parameter. 
The eigenvalues of the hopping matrix Q were determined by general purpose 
routines for non-hermitian matrices. Since the configurations are compara- 
tively small, direct determination of all eigenvalues and eigenvectors was 
possible throughout. 
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0.1 


0.5 


1.0 


1.5 


2.0 


3.0 


4.0 


5.0 




0.332 


0.314 


0.296 


0.286 


0.276 


0.267 


0.262 


0.260 



Table 1: Values of k for different (3. 

2.3 Topological charge in the continuum and on the 
lattice 

In this subsection we discuss the topological charge in the continuum and its 
geometric definition on the lattice. We state the Atiyah-Singer Index Theo- 
rem and the less known Vanishing Theorem which holds for QED2. 

We start with collecting some results for QED2 in the continuum which 



can be found in the classical papers |22|. There space-time is compactified 



on a sphere by requiring the gauge fields to approach a pure gauge at infinity. 
The case of QED2 on a (continuous) torus is discussed in ||23|. In both cases 
one can define the topological charge (Pontryagin index) 

iy[A] = d\Fu{x) , (6) 

which is an integer valued functional for classical (differentiable) gauge field 
configurations. The Atiyah Singer Index Theorem relates the topological 
charge to the index of the Dirac operator, which is given by the difference 
of the number of positive (n+) and negative (n_) chirality zero modes. The 
zero modes are eigenstates of the Dirac operator 

i]3ip{x) = i'jf.ldf, - ieAf,{x)]ip{x) = Eip{x) , 

with eigenvalue E = 0. Since 73 anticommutes with the zero modes can 
be chosen as eigenstates ip+yi^- of 73 with 'J31P+ = +"0+ cind 73V'- = — V'-- 
n+ (n_) denotes the number of independent states ip+ The Atiyah 

Singer Index Theorem ^ then reads 

u[A] = n+ - n_ . (7) 

For the case of QED2 there holds another index theorem, which is sometimes 



referred to as Vanishing Theorem |22 



u[A] > =^ n_ = 
6 



u[A] < ^ n+ = . (8) 

The Vanishing Theorem thus states that for non-zero topological charge ei- 
ther only positive or only negative chirality eigenmodes exist. 

When computing the topological charge for the lattice gauge fields, we 
use the geometric definition which is based on Liischer's idea of associating 
a principal bundle to each lattice configuration and defining its topological 



charge through the topological charge of the bundle [0. The case of QED2 
was worked out in [^. One obtains 



Mu] = ^E^p(^) e (9) 

The plaquette angle 6p{x) is introduced as (compare (||)) 

ep{x) = Im \nUp{x) , (10) 

and restricted to the principal branch 9p{x) G (— vr, vr). Note that configura- 
tions where Up{x) = —1 for some x are so-called exceptional configurations 
and Liischer's definition does not assign a value I'llU] to them. (This nota- 
tion should not be confused with the widely used terminology where one calls 
exceptional configurations those gauge field configurations with bad conver- 
gence properties in the inversion algorithm for the Dirac operator.) However, 
those configurations are of measure zero in the path integral. 

3 Analytic results for the eigensystem of the 
fermion matrix 

This section is devoted to some analytic results on the general structure of 
the spectrum and the chiral properties of the eigenvectors of the fermion 
matrix. These results are valuable tools for our later numerical investigation 
of the spectrum and are necessary for a proper interpretation of an eventual 
approximation of the continuum index theorems by the lattice model. 
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3.1 The general structure of the spectrum 

In this subsection we collect analytic results {SI - S4) for the spectrum of the 
fermion matrix M which partly can be found in the literature ^ Here 
we carefully distinguish between even and odd values of L since the structure 
of the spectrum changes according to L even or odd. A proper investigation 
of this fact is necessary for understanding an eventual realization of the index 
theorems (|^, (||) on the lattice. 

We start with noting that the hopping matrix Q, defined in @ is similar 
to its hermitian adjoint 

Tagrs = , (11) 

where 

Ts = Isues ® 73 with = I , = Ts . (12) 

This implies that a is a simultaneous eigenvalue of (eigenvector v), and 
of Q (eigenvector T^v). Let A be some eigenvalue of Q; then from general 
theorems A is an eigenvalue of and thus also of Q. 

SI : Eigenvalues of Q and hence also of M (compare (Q) ) are either real or 
come in complex conjugate pairs. 

Furthermore for L even Q is also similar to — Q 

SgS = -Q , (13) 

where 

E{x,y) = i-lT^'^^^ 6^,,y,6^,^y, Ispinor with = 1 , = E . 

(14) 

For L odd (|13]) does not hold, since the crucial condition (compare (|^)) 

^x±fi,y(y 1)" ' y y ^x±fi,y ; 

is violated for = 1 (or L). The property = I implies that if A is 
an eigenvalue of Q with eigenvector v, then also —A is an eigenvalue of Q 
(eigenvector Ev). Thus the spectrum obeys: 

S2: For even L, eigenvalues of Q come in pairs A, — A. 
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Let us now turn to the eigenvalues /i of the fermion matrix M itself. From 
one finds: 

S3: The eigenvalues /i of the fermion matrix M are related to the eigenvalues 
A of the hopping matrix Q through ^ = 1 — kX . 

This has the direct consequence that for any positive, real A there is a kq 
such that on that gauge field configuration /u(ko) = 0. 

By establishing a bound on the norm of Q, further information on the 
spectrum can be obtained. With (|^), using the fact that l + 'ju and l — 'ju are 
proportional to orthogonal projectors and |?7j,(a;)| = 1 it is straightforward 
to show 

II n II II II o 

oo = sup— ^ r— = 2 . 

The norm || .. || is defined to be the norm obtained by summing over all 
lattice and spinor indices and g is some (spinor) test function on the lattice. 
This implies for an eigenvalue A of Q with eigenvector v 

II II 

l-^l — ~n ii~ —HQ lloo = II Ql + Q2 ||oo ^ II Ql ||oo + II Q2 lloo = 4 . 

II ^ II 

Together with S3 one finds: 

S4: The eigenvalues n of M are distributed inside a circle with radius = Ak, 
with center 1 in the complex plane. 

The results SI - S4 have important implications for an approximate realiza- 
tion of the index theorems (|^), (§) on the lattice. For example a configuration 
with |t'z[?7]| = 1 lets one expect (due to (0) and (^) exactly one approximate 
zero-mode, i.e. exactly one eigenvalue of M with small modulus. This eigen- 
value is then necessarily real, since if the lowest lying eigenvalue was complex, 
both this eigenvalue and its complex conjugate would have the same absolute 
value. The emergence of one small and real eigenvalue should be invariant 
under small deformations of the background field U. Finally we remark, that 
S4 implies that for k < 1/4 (which corresponds to bare mass m > 0) exactly 
vanishing eigenvalues are excluded, and real eigenvalues are positive. This 
corresponds to the fact that for non-vanishing quark masses m also in the 
continuum exact zero modes are excluded. 
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3.2 Chiral properties of the eigenvectors 

In this section we prove a theorem on the chiral properties of the eigenvectors 
of the fermion matrix. It establishes that indeed, as already conjectured in 
the end of the last paragraph, only the eigenvectors with real eigenvalues can 
be interpreted as approximate zero modes. 

In the continuum the zero modes can be chosen as chiral eigenstates, since 
73 anticommutes with the Dirac operator. On the lattice the Wilson-Dirac 
operator has terms that explicitly break the chiral invariance and it does 
not anticommute any longer with 73. However, it is interesting to study the 
pseudoscalar density matrix 

Xihj) = {vi,'r3Vj) = vlV^Vj , (15) 

where Vi, Vj are the right eigenvectors of the hopping matrix Q with eigen- 
values Aj, \j. Due to (0) Vi, Vj are also eigenvectors of the fermion matrix M, 
with eigenvalues fii, fij related to Xi, Xj through S3. Since r3 is hermitian, 
we find 

Xij,i) = Xihj) and xihi) ^ R, (16) 
which shows that x is a hermitian matrix. The following theorem holds: 

Theorem 1 (Vanishing entries of the chiral density matrix): 

A, ^ a; ^ xi^,J) = 0, (17) 
and in particular for the diagonal elements 

Xi ^ R ^ X{i,i)=0. (18) 

Note that using S3 the conditions Aj 7^ A^ and Aj ^ R can immediately be 
written into conditions [li ^ JT] and /ij ^ R for the eigenvalues /Xj of the 
fermion matrix M. 

Proof: We define the matrix H = T^Q. From (|11|) it follows that H is 
hermitian and thus {v, Hv) G R for arbitrary vectors v. It follows (make use 
of the eigenvalue equation Qvi = XiVi) 

R 3 ivi,Hvi) = {vi,TsQvi) = Xi {vi,T3Vi) = Xixihi) ■ 
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Since xih'^) G R we conclude that whenever xih'^) is non-vanishing, Aj has 
to be reaL This proves (|l^) for the diagonal elements. The result for the 
diagonal entries of the pseudoscalar density matrix may be found in PT]]. 
Next consider 

R 3 {vi + Vj,H[vi + Vj\) = Xixii,i) + + + XjXihj) , 

U 3 {vi + iVj,H[v^ + iVj]) = Xixihi) + XjXUJ) - iKxij,^) + ^^^jXihj)- 

Since we already showed that the first two terms on the right hand sides are 
real (either Aj is real or x{hi) vanishes) the last equation implies (use (|TB|)) 

XiXihi) + XjX{h3) e R, 
- i\ xihj) + iXjXihj) e R. (19) 

Introducing the abbreviations xih j) = + = Xj + iyi and Aj = Xj + iyj 
the vanishing of the imaginary part in (|T9|) can be expressed as 



' b a' 








0" 


—a b 




. Vi + Vj . 




.0. 



The matrix on the left has determinant a? + b'^ . A non-trivial solution for the 
vector {xi — Xj,yi + yj)'^ is possible only if this determinant vanishes. This 
implies that whenever + 6^ 7^ 0, we must have 

•^i -^j ) Hi Vj Xi Xj . 

The condition + 6^ 7^ is equivalent to x{h j) thus (0) is proven. 
□ 

Theorem 1 has important consequences for an eventual interpretation 
of eigenvectors of the fermion matrix as candidates for approximate zero 
modes. In particular the statement (^) on the diagonal entries x{h i) of 
the pseudoscalar density matrix shows that only eigenvectors Vi with real 
eigenvalues /ij allow for non-trivial chiral properties. If one assumes that 
the index theorem of the continuum can be recovered from the continuum 
limit on the lattice, eigenvectors Vcompiex with (truly) complex eigenvalues 
are already ruled out as candidates for zero modes, since they always obey 
{^complex, ^sVcompiex) = 0, whcrcas normalized zero modes in the continuum 
obey (^',73^') = ±1- 
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4 Numerical results for the eigensystem of 
the fermion matrix 



In Section 3 we explored some general aspects of the spectrum and the chiral 
properties of the eigenstates of the fermion matrix using analytic methods. 
Here we numerically analyze the behaviour of the spectrum in Monte Carlo 
generated background configurations. 

4.1 Spectrum in fixed baclcground configurations 

In this subsection we discuss the spectrum of the fermion matrix using sev- 
eral fixed background configurations. We formulate some conjectures on the 
general structure of the spectrum, which will be tested in subsequent sections 
using the whole Monte Carlo ensemble of configurations. 

When the couplings approach their critical values (compare Section 2.2) 
we find p8[ that the physical signature (spectrum of triplet and singlet cur- 
rents, chiral condensates) of the lattice model is in agreement with the con- 
tinuum results for QED2 with massless fermions. Thus we expect that the 
eigenvalues of the full fermion matrix (including the constant term) resembles 
the spectrum of the massless continuum Dirac operator. We remark that we 
restrict ourselves to lattices L x L with even L in order to be able to use the 
maximal symmetry including S2. 

In Figures |I|.a - |l].f we show the spectrum of the hopping matrix Q for a 
set of 'typical' background configurations with various values of z/; (compare 
also |T0| , |11|). The Monte Carlo generated configurations are in equilibrium 
for the corresponding values of k and f5. The values of k were chosen near 
the critical k for the corresponding f3. 

It is obvious that all spectra obey the symmetry properties SI, S2, S4. 
The eigenvalues are either real or come in complex conjugate pairs, they are 
enclosed in a circle of radius 4 around the origin in the complex plane. Their 
distribution is symmetric with respect to reflection at the imaginary axis, as 
expected from S2 for even L. 

We remark, that due to S3 the largest real eigenvalues of Q are to be- 
come the smallest real eigenvalues of M and thus are the candidates for the 
approximate zero modes. The spectra for topological charge vi = 1 (Figs. 
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|l|.b and |l|.e) show one large real eigenvalue of Q (corresponding to one small 
real eigenvalue of M) as was predicted already in the discussion at the end 
of Section 3.1. This eigenvalue has the clear interpretation of corresponding 
to an approximate zero mode. 

In the sector with z/; = 2 (Fig. |l].c) the situation is more involved. The 
spacing of the two largest real eigenvalues is of similar size as their dis- 
tance from the nearest pair of complex eigenvalues. Thus the size of the 
corresponding eigenvalues of M is not a proper criterion to identify them as 
approximate zero modes. However from the analysis of the chiral properties 
of the eigenvectors in Theorem 1 it is clear that the 2 small real eigenvalues 
of M correspond to the approximate zero modes. 

It is interesting to note, that in all plots the number of real eigenvalues 
is equal to 4 | z/; | . It is generally believed ||^, |1^, |ll| that the smaller real 



eigenvalues are the 'would be' zero modes from the doublers, which are shifted 
to higher values by the Wilson term. Thus one is tempted to conjecture (for 
L even) 

# of real eigenvalues = 4|z/i| . (20) 



However, it has to be remarked, that (|20D cannot be exact for arbitrary gauge 
field configurations. A counter example is given by 

U2{xi,L) = —l Vxi = l,...L ; all other [/^(xi, X2) = 1 . (21) 

This configuration has Up{x) = 1 for all x G A and thus i^ilU] =0. It 
furthermore reduces the fermion matrix to the case of free fermions with 
periodic boundary conditions. The spectrum for this case can be computed 
using Fourier transform and is given by (ni,n2 = 1, ... L) 



1 _ 2^ l^cos (^) + COS (^) ± y 2 - cos2 (^) - cos2 

It is obvious that the imaginary parts vanish whenever ni,n2 G {L/2,L}. 
This establishes the fact that there are real eigenvalues although I'llU] = 0, 
implying that (pO|) can be violated for certain configurations. This possible 
violation was also confirmed, when we performed a numerical analysis of the 
spectrum for randomly chosen background configurations. In this study with 
static background configurations it also turned out, that the smoother the 
gauge field, i.e. plaquette variables close to 1, was chosen, the less likely a 
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Figure 1: Spectrum of Q in the complex plane for a 8 x 8 lattice. The 
background configurations in (a-c) were taken at f3 = 2.0 and n = 0.276. Plots 
(a), (b) and (c) correspond to i^i = 0, 1 and 2, respectively. The parameters 
for (d) and (e) are P = 4.0, k — 0.262, with (d) and (e) showing 1^1 = and 1. 
Finally in (f) we plot the spectrum, for [/ — 1^ k — 0.25 which corresponds to 
free massless fermions (note that we use mixed periodic boundary conditions 
and that the eigenvalues in (f) are degenerate). 
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violation of (^) became (although isolated configurations such as ( pT]) with 
all plaquette variables equal 1 still violate (^)). We conjecture that up to 
isolated configurations, is indeed correct for sufficiently smooth gauge 
field configurations (compare |10], 0). 



4.2 The relation between the topological charge and 
the real eigenvalues 



In the last section it was conjectured, that ( [20| ) should become correct for 
gauge field configurations which are smooth in some sense. In particular one 
is interested in the behaviour at large /3 since one wants to construct the 
continuum limit when /3 — >■ oo. It is known that with increasing (3 the gauge 
field becomes smoother (plaquette variables closer to 1) since local fiuctua- 
tions are stronger suppressed by the action. 

In order to test the conjecture that the probability p(/3) of finding (pO|) 
correct increases with /?, we plot p{P) in Fig. |^ as a function of P for various 
lattice sizes (L = 4,8,12, and 16). For each f3 the value of n was always taken 
to be approximately Hcrit for this P and lattice size, as discussed in Sec. 2.2. 
The values of /? are 0.1, 0.5, 1.0, 1.5, 2.0, 3.0, 4.0, 5.0. For L = 4 and 8 we 
used 10^ configurations and 10^ for L = 12 and 16. 

The figure shows clearly, that already for modestly high values of (3, the 
configurations obeying (pO| ) entirely dominate the path integral. For jS > 3, 
p{P) is essentially equal to 1 for all lattice sizes we analyzed. We remark, 
that we found only a weak dependence of p(/?) on k. 

What is interesting to note is the difference between the various lattice 
sizes. For fixed (3 < 3 the larger lattices lag behind smaller ones in obeying 



(20). This poses an interesting question: Do the curves for p{f3) reach a 
limiting curve when L increases? Ideally when performing the continuum 
limit one would first like to send L ^ oo followed by the limit (3 ^ oo. If 
a limiting curve with p{f3) = 1 for all (3 larger than some finite bound /3o 
emerges one could work with (^) for all P > Pq, independent of the lattice 
size. If no such limiting curve emerges, the bound Po{L) depends on L. 
However, for fixed L we can use (pO|) for sufficiently large p. 

There is a second detail which deserves discussion, but is peculiar for 
QED2. Although the classical index theorems (0), (^) are used only as a 
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G O 4x4, 10 confs. 

Q □ 8x8, lO'* confs. 

O—o 12x12, 10^^ confs. 
A-----A 16x16, 10^ confs. 
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Figure 2: Probability of finding ( [2D| ) correct as a function of p. We show 
our results for lattice sizes L = 4,8, 12 and 16. The symbols are connected 
to guide the eye. 



guideline in our decomposition of the path integral, the distinct role of the 
trivial sector in the Vanishing Theorem (which only holds for QED2) 
should be further analyzed. The theorem makes no statement on the trivial 
sector and only the Atiyah Singer Index Theorem (0) applies. Thus there is 
no immediate hint, that (^Oj) should also hold for ui = 0. In order to settle 
this issue, we show in Fig. |] our results for decomposed with respect to 
1 1// 1 on a L = 16 lattice. The raw data for this plot are the same as in Fig. ^ 
We find, that for /? > 1 the different sectors are essentially uniform in 
obeying (|20|) . In particular we find no anomaly for the trivial sector. We 
remark, that for larger f3 and higher |z/i| it is difficult to obtain good statistics, 
since the higher sectors already have a considerably smaller weight in the 
path integral (compare Section 5.1). The observed high values of p(/3) for 
/3 < 1, > 1 might be explained by the increase of the effective action with 
|z//| (see Section 5.1). Increased values of the effective action for higher |z//| 
damp additional quantum fluctuations of the gauge field in the higher sectors 
leading to smoother configurations with enhanced p(/3). When (3 increases, 
the configurations are forced to be smoother also in the lower sectors leading 
to the observed uniform behaviour in \i>i\ for higher f3. 
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P 

Figure 3: p{P) decomposed with respect to |z/;|. The data were computed 
from the configurations on a 16 x 16 lattice aheady used in Fig. |^. The 
symbols are connected to guide the eye. 

4.3 Properties of the spectrum in ensembles of back- 
ground configurations 

Pursuing the program outlined in the introduction one would like to go 
over from individual background configurations to an investigation of the 
behaviour of the spectrum in a whole ensemble of configurations with fixed 
winding number. In this section we establish some properties of the average 
behaviour of the spectrum in the ensemble of the path integral. In particular 
we concentrate on the real eigenvalues which are intimately connected to the 
topological charge, as established in the last section. 

We start with showing a schematic picture (Fig. ^) for the distribution 
of the real eigenvalues of the hopping matrix Q. 

The pattern displayed in Fig. ^ emerges clearly for smooth configurations 
and becomes more pronounced for increasing values of (3 (see also [0). We 
find that the real eigenvalues are concentrated in the small domains A,B,A* 
and B*. The eigenvalues in A are related to the eigenvalues in A* through 
the symmetry A ^ —A (S2) and the same holds for B and B*. We find that 
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Figure 4: Scliematic picture of tlie distribution of the real eigenvalues of the 
hopping matrix Q. The shaded areas are the domains A*,B*,B and A where 
the real eigenvalues concentrate. In the top line we indicate how the sign of 
the matrix elements xih'^) corresponding to the eigenvalues in the domains 
is related to the sign of the topological charge vi. 



the corresponding diagonal entries of the pseudoscalar density matrix xih'^) 
are equal (as a matter of fact this can be proven analytically, using similar 
techniques as in the proof of Theorem 1). Also the sign of the diagonal 
entries obeys a simple pattern: It coincides with the sign of z// for the regions 
A, A* and equals minus the sign of vi for B and B* . This shows that the 
eigenvectors corresponding to domain A, which are the approximate zero 
modes (their eigenvalues of M approach 0) also have the correct chirality as 
expected from (|^), (§). 

We find that with increasing (3 the domains A and B shrink and are 
shifted towards the limiting values 4 and respectively (see also JTO], In 
order to quantify this statement, we computed the average values A, B and 
the standard deviations aA-, c"b of the real eigenvalues in each of the domains 
A and B. The eigenvalues can be attributed to A or i? due to the chiral 
properties of the eigenvectors (i.e. the sign of x(i,2)). 

Fig. ^.a shows our results for 8 x 8 lattices and various values of (3. For 
/3 < 2 we evaluated 10^ configurations in non-trivial sectors, for /3 = 3, 4, 5 
we used 10^ nontrivial configurations. The symbols give the average values 
A (diamonds) and B (squares) of the real eigenvalues in A,B. horizontal 
bars show the corresponding standard deviations a a, ctb- Note that this is 
not some kind of errorbar giving the size of statistical errors, but gives the 
average size of the regions A and B. It is obvious, that with increasing f3 the 
real eigenvalues get confined in smaller regions closer to the limits and 4. 

The same picture also holds for each sector separately. In Fig. |^.b we show 
the same quantities as in Fig. |^.a, but restrict the ensemble to configurations 
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Figure 5: Average behaviour of the distribution of the real eigenvalues for 
various values of (3. The data were taken on a 8 x 8 lattice (see the text for 
statistics). Fig. ^a (left-hand plot) shows the average values A, B of the real 
eigenvalues in the domains A (diamonds, right branch) and B (squares, left 
branch). The horizontal bars show the size of the standard deviations a a-, c^b 
from the average values. We stress that this is not a statistical error, but 
gives the average size of the domains A and B. In Fig. ^.b (right-hand plot) 
we show the same information but take into account only configurations with 
1^1 = 1. 
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with 1/1 = 1. Again we find that with increasing f3 the real eigenvalues get 
confined in smaller regions closer to the limits and 4. 



^ 1 

^ 1 

^ 1 

^ 1 

-<) 1 

-o 1 

0.0 1.0 2.0 3.0 4.0 

real eigenvalues 

Figure 6: Average behaviour of the distribution of the real eigenvalues for 
different ui. The data were computed from 2300 configurations in non-trivial 
sectors on a 16 x 16 lattice at /3 = 2, k = 0.276. The symbols give the average 
values A, B of the real eigenvalues in A (diamonds, right-hand branch) and 
B (squares, left-hand branch). The horizontal bars are the corresponding 
standard deviations aA^crs- 

If one compares the different topological sectors (/? fixed) one finds that 
the average values A, B of the real eigenvalues in A, B are rather indepen- 
dent of vi and also the size of the domains A, B is rather invariant. This is 
demonstrated in Fig. ^ where we show average value and standard deviation 
of the eigenvalues in A and 5 as a function of the topological charge. We 
used 2300 configurations in non-trivial sectors, generated on a 16 x 16 lattice 
at /5 = 2,fi: = 0.276 (^ KcHt)- 

Also the diagonal entries of the pseudoscalar density matrix show 

a simple behaviour, monotone in [3. We computed the average values Xa-,Xb 
of the x{h cind their standard deviation o"^ from this value for each domain 
A, B separately. Fig. |^ shows our result for the same sample used in Fig. ^.b. 
z// is positive (=1), and as expected from (0) and (|) domain A (diamonds) 
has positive xi^i i)i while in region B (squares) we find negative values. Both 
branches show a clear trend towards the continuum values —1 and 1 with 
increasing (3. Also the interval containing xih'^) (horizontal bars) shrinks for 
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Figure 7: Average values Xa,Xb of the entries xI^jO of the pseudoscalar 
density matrix corresponding to eigenvalues in A (diamonds, right branch) 
and B (squares, left branch) in the z/; = 1 sector. The horizontal bars show 
the size of the standard deviation from the average value. The data were 
computed from the same configurations used in Fig. ^. 

larger /3. 

We summarize the emerging picture for the structure of the spectrum in 
the ensemble of the path integral as follows {L is even): 

El: Only real eigenvalues can correspond to approximate zero modes, since 
for complex eigenvalues the diagonal entries of the pseudoscalar density ma- 
trix vanish identically (Theorem 1). 

E2: For large (3 the path integral is entirely dominated by configurations 
obeying # of real eigenvalues = 4|i/;| (Fig. 

E3: The real eigenvalues group according to the scheme depicted in Fig. ^. 
The domains A, B, A*, B* are related through S2 and the doubler symmetry. 

E4: With increasing (3 the domains A and B shrink and move towards the 
limiting values 4 and (Fig. ^). 

E5: The sign of the diagonal entries x{h i) of the pseudoscalar density matrix 
coincides with the sign of vi for the eigenvalues in A, A* (as expected from 
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(0), (^)) and is minus this sign for B,B*. With increasing [3, the values of 
xih'^) g^'t confined in shrinking regions approaching the limiting values ±1 
(Fig. 0). 

E6: The behaviour of the eigensystem is essentially uniform in ui (Figures ^ 
and 1^. 

In the remaining sections we use these results for analyzing the contri- 
bution of different topological sectors to the path integral and to vacuum 
expectation values of various operators. 



5 Applications 

In this section we bring in the harvest from the previous investigation of the 
interplay between topological charge and behaviour of the spectrum of the 
fermion matrix. We will demonstrate, that the analytic results SI - S4 and 
Theorem 1, together with the results El - E6 on the spectrum in ensembles 
of background configurations provide a powerful tool for analyzing physical 
questions in lattice gauge theory. Here we concentrate on local and bulk 
quantities. The study of the dependence of the spectrum (which requires 
the computationally more demanding analysis of two-point functions (see 
also [|^)) on the topological charge is reserved for a forthcoming publication 



5.1 Effective action and fermion determinant 

In this first application we concentrate on the behaviour of the effective ac- 
tion and in particular of the fermion determinant in the various sectors. 

After integrating out the fermions (note that we use 2 explicit flavors of 
fermions (not to be mixed up with the doublers implicit in the Wilson action) 
one obtains the effective action for the gauge fields 

S'ff[U] = Sg[U] + Sf^[U] = Sg[U] - lndetM[t/]2. (22) 

It is interesting to analyze typical values of the effective action in various 



topological sectors. For the gauge field part of the action it is known 



that there exist lower bounds of the action for each sector. In a previous 
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article |30[ we analyzed the role of these bounds in the case of pure U(l)2 
gauge theory on the lattice (we will comment on this below). Due to the 
relation of topological charge and zero modes via the index theorems, also 
the fermion determinant depends on the topological sectors in a non-trivial 
way. 

In Fig. P we show our results for the averages of the gauge field action 
Sg{i') and the effective action from the fermions S ^ [v) in each of the sectors. 
We give the results per plaquette, i.e. normalized with L~^. 
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Figure 8: Average value of the contributions of the gauge field action 
Sg{vi)/ L"^ ((a) left-hand side) and of the fermions S^j \vi) j L'^ ((b) right-hand 
side) to the effective action. The average was taken for each topological sec- 
tor separately. The data were obtained from a simulation on a 16 x 16 lattice 
at /3 = 2, K = 0.276 using 10^ configurations. In (a) we connect the symbols 
to guide the eye. The straight line in (b) was drawn using ( P^ ) below. 

Fig. |[a (left-hand plot) shows the average value Sg{y) of the gauge field 
action in each topological sector separately. We find a behaviour which is 
essentially linear in |z/;|. If one compares this result with the pure gauge 
theory |30| the linear behaviour is rather surprising. In the quenched case 
we found a behaviour of the gauge field action which is quadratic in vi with 
a variation in vi which is of one order of magnitude smaller than here. This 
result was understood in a quantitative way using an argument based on the 



topology of the configurations [|30[. The essentially linear behaviour (plus 



eventually a small quadratic term) in Fig. ^.a has no such simple explanation. 
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In the unquenched model it is not the gauge field action alone which is 
minimized, but the effective action (^2]) which is the sum of the gauge field 
action and the contribution from the fermion determinant. Obviously the 
influence of the fermions adds a linear term which entirely changes the picture 
of the pure gauge theory. The interplay of the two terms in (p2D does not 



allow for a simple topological explanation as in |^ 

There is however a lesson from the study of the pure gauge theory pO 
which is also conflrmed in the unquenched model: The lower bounds for the 
action in each topological sector are at least one order of magnitude smaller 
than the average value of the gauge fleld action. We also computed the dis- 
tribution of the values of the action, and flnd that quantum fluctuations keep 
most of the conflgurations high above the lower bounds in each sector [pp]]. 
The same observation is true for the model with fermions. This is in particu- 
lar remarkable, since the lower bounds given in [^] are saturated by classical 
conflgurations with constant electric fleld, carrying topological charge. The 
path integral for the Schwinger model on the torus including all topologi- 
cal sectors can be constructed by dressing these classical conflgurations with 



quantum fluctuations |23|. The observed high values of the gauge fleld action 
for the main part of the conflgurations is a demonstration of the importance 
of quantum fluctuations (see for a more detailed discussion). 



Fig. |[b (right-hand plot) shows the average value S^^ {yi)lL'^ of the 
contribution of the fermions to the effective action for each topological sector. 
Again we flnd a behaviour essentially linear in |z/;|. This behaviour can be 
understood in a quantitative way using the results from the previous sections, 
as will be shown now. 

From S2, S3 it follows [L even) that real eigenvalues p of M come in pairs 

p(l) = 1-KW , p(^) = l + KW , (23) 

with < w < 4. E2 and E3 establish, that for large enough (3 

of real pairs = 2|z/;| , (24) 

and that the real pairs are equally distributed among the domains AU A* 
and BUB*. The determinant of M is the product of all eigenvalues. For the 
contribution of the real pair (ESI) one obtains 



R = p(^)p(^) = 1 - (^wf. (25) 
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In E4 we establish, that the regions A and B (and of course also their S2 
images A* and B*) shrink and approach their limiting values and ±4). Thus 
we obtain characteristic contributions from the real eigenvalues in AUA* and 
B U B* described by the average values A, B 

Ra = I- {kAY , Rb = 1- {kBY . (26) 

The contribution of the complex eigenvalues is more difficult to handle. How- 
ever by inspecting the spectra for several configurations (see Fig. |l| for a few 
examples) we find that whenever increases by one unit, one of the com- 
plex quadruplets vanishes and is turned into two real pairs. The rest of the 
complex eigenvalues remains rather unperturbed giving some positive (note 
that complex eigenvalues come in complex conjugate pairs) constant c. As 
long as \i>i \ is considerably smaller than the overall amount of eigenvalues 2L^, 
we make the following ansatz for the contribution of the fermion determinant 
to the effective action: 

Sff = -lndet(M)2 ~ -ln(c(i?Ai?B)''''')^ = C - 2\n\RARB\ ■ 

(27) 

The expression gives a behaviour linear in and predicts also the slope of 
the line. The constant term has to be taken from the numerical data. Eval- 
uating 2300 configurations for these parameters {L = IQ, (3 = 2, n = 0.276) 
we obtained A = 2.996(8), B = 0.421(7). The straight line in Fig. |.b was 
drawn using ( pT]) with these numbers. We find good agreement with the 
Monte Carlo data. Eq. (^) slightly overestimates the Monte Carlo results 
for larger \i>i\. This might be due to the fact that by increasing \i>i\ by one, 
four new real eigenvalues emerge causing one complex quadruplet to vanish. 
Thus the contribution from the complex eigenvalues is not entirely indepen- 
dent of I 1 as we assumed in our ansatz. This effect should become smaller 
with increasing L. We believe that the essentially linear relation, and the 



good performance of (|2^) shows that the behaviour of the fermion determi- 
nant as a function of z// can be understood from the analysis of the spectrum 
of the fermion matrix in Sections 3 and 4 in a quantitative way (compare 
also iTTI). 



We end this section with a remark on the positivity of the fermion deter- 
minant det M. From SI - S3 we know that the complex eigenvalues come in 
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quadruplets 



(1) 



1 — K, {x + iy) 
1 + K {x + iy) 




1- K{x-iy) , 
1 + K {x — iy) , 



(3) 



/i 



contributing a factor 



4 



(1 - Kxf + (1 + fixf + K^y"^ 



to the fermion determinant. Obviously the factor from a complex quadruplet 
is strictly positive. 

The second order polynomial 1 — {nwY which one obtains for the contri- 
bution of a pair of real eigenvalues (|25|) becomes negative for kw > 1. Thus 
a negative value of the fermion determinant can only come from pairs of real 
eigenvalues. Since the emergence of real eigenvalues is connected to gauge 
field configurations with non vanishing topological charge (compare E2) the 
positivity properties of the determinant are as well. In the literature one 
often encounters the approximation 



used to simulate odd numbers of flavors. This might possibly lead to a wrong 
weight from non-trivial topological sectors. 

In order to investigate this possibility we performed several simulations 
on L = 8 lattices at /5 = 2. The values of k where chosen above the critical 
ft (which is Kcrit ~ 0.276 for /3 = 2) in order to force the determinant to 
negative values, n varied between 0.280 and 0.335. For k > 0.320 we found 
negative values of the determinant in the non-trivial sectors. It is interesting 
to remark, that we then also found a considerably enhanced abundance of 
configurations in the higher sectors. The mechanism for this trapping is due 
to the contribution of the real eigenvalues to the effective fermion action, 
which is essentially the sum of the logarithms of the eigenvalues squared 
(note that we squared the determinant for the updating). Once one of the 
real eigenvalues assumes a negative value (thus causing the determinant to 
be negative), it can continuously go back to a positive value only by passing 
through zero. This however drives the logarithm and thus the effective ac- 
tion to infinity. The only way back is through a discontinuous jump, which 



det(M) 



det(M)| , 
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Figure 9: Time series of the fermion determinant and the topological charge. 



probably costs a big amount of action. This creates a barrier, which traps 
the system in a topologically non-trivial sector. 

Fig. y shows a time series for the determinant and the topological charge 
from 10^ configurations at (3 = 2, k = 0.33. For approximately the first 1000 
configurations the topological charge is mainly zero with a few configurations 
that fluctuate to ui = ±1. The determinant has only positive values. Then 
for approximately 3400 configurations, the determinant has negative values, 
and the topological charge is mainly trapped at ui = —1 implying that the 
fermion matrix has 4 real eigenvalues (and thus can have negative values). 
There are a few configurations with ui = or 2. From the sign of the 
determinant we conclude that there still are quadruplets of real eigenvalues, 
although e.g. ui vanishes. These are just the configurations violating (|20| ) in 
agreement with our observation in Fig. However, the system is obviously 
trapped in the ui = —1 sector for more than 3000 configurations. After that 
a tunneling back to positive values of the determinant occurs and the system 
evolves again with positive values of the determinant and the distribution of 
the values of ui is again symmetric around 0. 
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5.2 The pseudoscalar density 

In this section we discuss the contributions of different topological sectors 
to the vacuum expectation value of the pseudoscalar density. This operator 
plays an important role in the so called 'fermionic definition' |]^ - of the 
topological charge which uses the chiral Ward identity to define the topolog- 
ical charge and is also used in the lattice derivation of the Witten-Veneziano 



formula [31 



Due to translation invariance, the vacuum expectation value of the pseu- 
doscalar density can be written as 

(V^(x)73^(x)) = -l^(^(a;)73^(x)) = J— J [dU]e-''"^^^ Tr (M-'Ts) . 

(28) 

In Fig. |10|we show our results for the average value Tt{M~^T3) in ensembles 
of background configurations with different {L = 16, (3 = 2, k = 0.276, 10^ 
configurations). Note that in these averages, the observable was not multi- 
plied with the Boltzmann factor exp{—S'^^^U]). 

The Monte Carlo data show a linear behaviour in z// with a positive slope. 
Such a behaviour can already be expected from the spectral decomposition 
of the continuum Dirac operator in classical background configurations. The 
continuum operator i]3 is hermitian, implying that eigenstates to different 
eigenvalues are orthogonal. If ip is an eigenstate with eigenvalue E, then 
due to {i]3,73} = 0, 'j^ip has eigenvalue —E. Thus the matrix elements 
< ipl'y^l'ip > vanish whenever ip is not a zero mode. The only states that 
can contribute in the spectral decomposition of the pseudoscalar density are 
the zero modes. Since their number is equal to u due to Index and Vanish- 
ing Theorem (|^), the linear behaviour follows for classical background 
configurations. 

The lattice Wilson-Dirac operator is not hermitian, even worse, not nor- 
mal. Thus its (right) eigenvectors do not form an orthonormal basis, and 
the continuum argument cannot be taken over to the lattice. However, we 
established {El) that the real eigenvalues on the lattice take over the role of 
the zero eigenvalues in the continuum. We furthermore demonstrated {E2) 
that the number of real eigenvalues is proportional to |z/;|. Finally we found 
(E5) that the chiral properties (i.e. the sign of the entry in the pseudoscalar 
density matrix) of the real eigenvalues in the physical region A behaves as 
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Figure 10: Average values of Tr(M"^r3) in different topological sectors. 
The symbols show the results from the Monte Carlo simulation {L = 16, f3 = 
2,K. = 0.276,10"^ configurations). The symbols are connected to guide the 
eye. 

in the continuum. These three ingredients allow to understand the observed 
linear behaviour of the pseudoscalar density. 



6 Discussion 

In this study we have demonstrated, that the existence of an integer valued 
lattice equivalent of the topological charge has important consequences. We 
believe, that the use of the topological charge on the lattice goes far be- 
yond its applications in continuum formulations of gauge theories. We have 
shown, that a careful analytic investigation of the spectral decomposition of 
the fermion matrix supports the interpretation of the real eigenvalues (and 
the corresponding eigenvectors) as the zero eigenvalues (zero modes) of the 
continuum Dirac operator. It should be stressed again, that this is a result 
for all configurations contributing in the continuum limit, while the Index 
Theorems in the continuum only hold for classical gauge fields. When con- 
centrating only on the real eigenvalues (which is much less demanding than 
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treating the whole spectrum), numerical methods were successfully used to 
analyze their distribution when approaching the continuum limit. We have 
shown, that the results for the spectrum can be used to understand the in- 
terplay between topological charge and physical quantities in a quantitative 
way. 

As discussed in the introduction, QED2 has many features in common 
with QCD4. The physical picture we obtained for the 2-dimensional theory 
thus provides an interesting model for the 4-dimensional case. We remark, 
that not only the physical picture, but also some of the mathematical results 
we obtained may be extended to QCD4. In particular Theorem 1 which is 
the conceptual backbone of this study has a straightforward generalization 
to the 4-dimensional case 



Other aspects are certainly considerably simpler in the 2-dimensional 
world. In 4-dimensions and for non-abelian gauge groups, the fermion ma- 
trix has much larger dimension, and explicit diagonalization becomes a chal- 



lenging enterprise |22|- However, with the identification of the real part of 
the spectrum as the trace of non-trivial topological charge, related physi- 
cal questions might be tackled without knowing all of the spectrum. When 
concentrating only on the real part, computational simplifications might be 
possible [0 . 



Although from a conceptual point of view most satisfactory, the geometric 
definition of the topological charge is very costly to evaluate for non-abelian 
gauge groups. Other definitions such as the 'fermionic approach' using the 
chiral Ward identity 0-0 might be simpler to compute. The latter ap- 
proach essentially makes use of the pseudoscalar density operator. We expect 
that also the results for the pseudoscalar density in Section 5.2 can be taken 



over to QCD4 [12 1 and might lead to a considerable simplification of the 



fermionic ansatz for the topological charge. 
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